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PACS 12.40.Ee - Statistical models of strong interactions 
PACS 13.88.+e - Polarization in particle interactions 
PACS 25 . 70 . Mn - Fragmentation (nuclear reactions) 

Abstract. - With the parametrization of parton distribution functions (PDFs) of the proton 
by Softer et al, we extend the valence quark contents to other octet baryons by utilizing SU(3) 
flavor symmetry. We find the method practically useful. Fragmentation functions (FFs) are 
further obtained through the phenomenological Gribov-Lipatov relation at the x — > 1 region. Our 
results are compared with different models, and these different predictions can be discriminated 
by upcoming experiments. 
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i" 1 1 . Introduction. — A major property of the theory 
of strong interactions, i.e. quantum chromodynamics 
(QCD), is color confinement, i.e., quarks and gluons can- 

1^ not be free but must be bounded inside hadrons. Due to 
the complicated non-perturbative effects of strong interac- 
fion, nowadays, it is still difficult to calculate baryon struc- 

C^ture functions, in terms of parton distribution functions 
(PDFs), from the basic principle of QCD. One possible di- 
section is lattice QCD, however, its calculation is largely 
constrained by the computational capability. Therefore, 

T 7^ phenomenological models based on the spirit of QCD still 
play an important role in the studies of baryon structures. 



Statistical models, providing intuitive appeal and phys- 



x, 

^ ical simplicity, have made amazing success. Angelini and 
Pazzi first verified that valence quark distribution has a 
statistical behavior when the Bjorken variable x is larger 
than 0.1. Utilizing a thermodynamical form as an input to 
the first order QCD Q 2 -evolution, they obtained a good 
fit at 5 < Q 2 < 100 GeV 2 pQ. Cleymans et al. then 
argued that Pauli principle should have effects on PDFs, 
and based on the MIT-bag model they got a statistical way 
to generate reasonable PDFs [2]. Along this way, several 
other models were proposed [51-112] . where quarks follow 
Fermi-Dirac distribution and gluons obey Bosc-Einstein 
statistics. For example, Mac and Ugaz [3] incorporated 
first order QCD correction and Devanathan et al. [5] ex- 
hibited the scaling behavior. Afterwards, Bhalerao et 
al. [7j considered finite-size correction, and Trevisan et 
al. [Sj introduced a confining potential and took gluon 



splitting into account. Recently, Zhang et al. established 
a new and simple statistical model in terms of light-front 
kinematic variables and calculated the nucleon PDFs [T2] . 
Besides, models based on other statistical principles were 
also studied. For instance, Zhang et al. [13] constructed 
a model using the principle of detailed balance, and the 
Gottfried sum they got is in surprising agreement with ex- 
periments. Singh and Upadhyay |14| extended this model 
to have spin considered with modifications. Further, Al- 
berg and Henley [TS] tracked the detailed balance model, 
and obtained PDFs for the proton as well as the pion. 

BBS model, established by Bourrely, Buccella and Sof- 
ter [HHH], allows PDFs to be directly expressed in terms 
of x. By using 8 free parameters, the model can obtain 
all polarized and unpolarizcd PDFs of the proton, and it 
well fits to experimental data in a broad range of x and 
Q 2 . Extension of the model to other baryons is still in 
progress, but one of the limitations is in the determina- 
tion of unknown parameters due to the currently relatively 
lack of relevant experiments. We stress that, SU(3) flavor 
symmetry [16] presents a kind of connection of quark dis- 
tributions between octet baryons, and by using it, we get 
PDFs of all octet baryons from the BBS PDFs of the pro- 
ton. 

On the other hand, fragmentation function (FF), which 
describes the process of quark hadronization, is another 
basic quantity in QCD. One application of PDFs is to 
study the fragmentation function D q (z). At present, it will 
be very useful if there exists simple connections between 
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them, so that one can predict the poorly known D q (z) 
from q(x). One tentative attempt is the so-called Gribov- 
Lipatov relation [T7], where 



xq h {x) = D h q {z), 



(1) 



is suggested as a phenomenological approximate when z{= 
x) — Y 1. Utilizing this relation, we get FF ratios, from our 
calculated PDF ratios, u(x)/d(x) and u(x)/s(x), for p and 
A, at the z — > 1 region. Recently, Soffer et al. gave FFs 
out of the same statistical spirit, and they parameterized 
the free parameters for p and A directly from experimental 
data [T5j. On the contrary, our results are based on SU(3) 
flavor symmetry. 

The paper is organized as follows. First we briefly 
review the quark-diquark model, perturbativc QCD 
(pQCD), and BBS model. Then we extend the BBS PDFs 
to octet baryons through SU(3) flavor symmetry, and com- 
pare results among different models. 

Quark-diquark model. — SU(6) quark-diquark 
model, based on the SU(6) quark wave-functions, was ap- 
plied to nucleon structures when x is large, first by Feyn- 
man 19J. In this model, when a valence quark of baryon is 
probed, we can reorganize the other two quarks as a mix- 
ture of a scalar diquark with spin and an axial vector 
diquark with spin 1, denoted as S and V respectively. 

Therefore, the unpolarized PDF is written as 



q{x) 



c q a< 



(x) 



,(x), 



(2) 



where c q and c q , the weight coefficients, can be de- 
termined from SU(6) quark-diquark wave-functions, and 
their values are different for different baryons; a D (x) (D = 
S, V) represents the probability of finding the diquark in 
the state D, and when expressed in terms of light-cone 
wave- function [5Uj, is 



a D (x) (x [d 2 k±]\(p(x,k±)\'' 



(3) 



with normalization to 3 when integrated in x € [0, 1]. 

To calculate explicitly, we can use the Brodsky- 
Huang-Lepage (BHL) prescription for the light-cone wave- 
function in momentum space, 



(p(x, k±) cx exp 




, (4) 



where m q and m D are masses of quarks and diquarks, 
respectively, and the parameter a D ~ 330 MeV. 
Similarly, for a polarized baryon, the PDF reads 

Aq(x) = c q a s (x) + c q a v (x), (5) 
a D (x) = /[d 2 fc±]Mx,fc±)| 2 |W D (x,fc±)|, (6) 



W D {x,k±) = 



(fc++m,) 2 -fc 2 
(fc++m 9 ) 2 + fc 2 



(7) 



where k + = x^J (m 2 + k\)/x + (m 2 + fc 2 )/(! — x). 

The PDFs of five baryons calculated from the quark- 
diquark model are listed in Table [T] [3T] . 

Perturbative QCD . — In pQCD, the quark distri- 
bution in the hadron h satisfies the counting rule [22] . 



q h {x) cx {l-xf, p = 2n-l + 2AS z , 



(8) 



where n is the minimal number of spectator quarks, and 
AS Z = \S1 — equals and 1 for parallel and anti- 
parallel quark helicity with respect to the hadron [23] . 
Keeping to the next-to-leading order, 



B 3 X 



(1 



B, 



a 



« 3.-1/2 



Ba 



« 3-1/2 



(l-x)\ (9) 



i>6 



-1/2, 



1 



(10) 



where B n = B(l/2, n+ 1) is f3- function. For each baryon, 
there are 5 constraints among 8 parameters [21] , so only 
3 are left free. Parameters for all octet baryons are given 
in Table H 

A more brief form of pQCD (called brief pQCD here- 
after) , considering only the leading terms of quark helicity 
distribution, is expressed as [24) 



17 \ A H - 

1i 0) = -5- x 
B3 



1/2 (1-*) 3 , 



qf(x) 



G H „-l/2 

B, X 



(1 



(11) 



112) 



For the nucleon, 4 constraints are given. They are (1) nor- 
malizations for qi, i.e., Ni = A qi + C qi ; and (2) the corre- 
sponding polarizations, i.e., Aqi = A qi — C qi . Parameters 
can be extracted through S = Au + Ad + As « 0.3, and 
the Bjorken sum rule T p — T n = \qaI 9v ~ 0.2, from polar- 
ized DIS experiments [25]. Besides, all octet baryons are 
related to each other by SU(3) flavor symmetry. So there 
is no free parameter in the brief pQCD. All parameters 
are listed in Table |3] 

BBS model. — In the BBS statistical approach, the 
nucleon is assumed to be a massless parton gas in equilib- 
rium at a temperature T. The parton distribution p(x), 
at an input energy scale Qq, can be written as 



p(x,Q 2 ) oc [e< x - x <»V*±i\- 1 , 



(13) 



where "+" is for fcrmions (quarks and antiquarks) and 
"—"for bosons (gluons). Xq p is the thermodynamical po- 
tential of parton p, and x is the universal temperature. 
PDFs for u, d quarks of helicity h, are [5] 



xq h (x,Q 2 Q ) 



xq h (x,Q 2 a ) 



Ax" 



e (z-Xo q )/* + 1 e x / s + 1 ' 



MXo q h ) 



Ax b 



(14) 



(15) 
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where we have g, = qf + q~ , Aq t = qf - q^ . From a 
next-to-leading-order fit of a selection of 233 data points at 
Qq = 4 GeV 2 , parameters for the proton were determined, 

x = 0.099, b = 0.40962, b = -0.25347, (16) 
A = 1.74938, A = 0.08318, A = 1.90801, 
X~ d = 0.30174, X+ = 0.22775, 
X+ u = 0.46128, X m = 0.29766. 

Considering the asymmetry, PDFs for s, s quarks [TT] 
are formulated as 



xs h (x 



In (l + e kx °*/> 



; (*-x \)/a + i i n ^ + e fcX + «/*) 
A K x l 



xs h (x,Q 2 ) 



e x/x _|_ J ' 

A x 2b ° ln (l + e-^Vj 

A s x 



(18) 



e x/x _|_ I ■ 

For the proton, the parameters are, 



X+ = 0.08101, X t 



0s 



0.2029, 



2.05305, A s = 0.05762. 



(19) 



From above expressions and by utilizing SU(3) flavor 
symmetry (Table S}, we get PDFs for all octet baryons 
as shown in Fig. 1 to Fig. 5. 

Fragmentation functions. — Fragmentation func- 
tion for u, d, s quarks, given by Soffer and Bourrely [18j . 
is of the form, 



D B q {x 7 Ql) 



A B q X*x b 
,(*-*?)/* + 1 : 



(20) 



where X^ is the potential related to the fragmentation 
process q — >■ B. They then determine the parameters 
as UM 



d = 0.428, b = 0.200, x = 0.099, (21) 



zlA _ /t A 
A u — A d 

Al = A p d = 0.264, A? = 1.168, Af = 1.094, 

XI = 0.648, Xf = 0.247, 

X^ = 0.296, X S A = 0.476, 

where the input energy scale was fixed at Qo = 0.632 GeV, 
and the yielded x 2 = 227.5 for 206 experimental points. 

Comparisons between models. We calculate 
PDFs as a function of x for models mentioned above, and 
our extension of BBS model based on SU(3) flavor sym- 
metry is compared with them. Comparisons are shown in 
Figs. 1,2,3,4,5. 



Below we discuss physical issues at the x — > 1 region. 

In Fig. 1, we illustrate Au v (x)/u v (x) of the proton. 
We can see that quark-diquark model, pQCD model, and 
brief pQCD model all give Au v (x)/u v (x) —> 1 when 
x — } 1, that is, the proton is completely positively po- 
larized. However, BBS model gives a different result of 
0.779. Fig. 1 also shows Ad v (x) /d v (x) , d v (x)/u v (x) for 
the proton. For Ad v (x)/d v (x), pQCD model and brief 
pQCD model give Ad v (x)/d v (x) — > 1 when x — > 1, while 
quark-diquark model and BBS model predict negative val- 
ues. For d v (x)/u v (x), the prediction of BBS model is 
close to those from pQCD and brief pQCD models, while 
the result of quark-diquark model is different from that of 
other models, i.e., it gives a larger value when x — > and 
approaches to when x — > 1. 

Fig. 2 gives the results of A. As is shown, for 
Au v (x)/u v (x), SU(3)-extended BBS model gives much 
smaller prediction than the other models, while for 
As v (x)/s v (x) and u v (x)/s v (x), predictions are close to 
each other. Such result means that our extension of BBS 
model through SU(3) flavor symmetry is reasonable in the 
studies of baryon structures. 

The comparisons for other baryons are also obtained. 
In Figs. 3 and 4, we make comparisons for E + and 
S°. We can see that in SU(3)-extended BBS model, 
Au v (x)/u v {x) approaches to 0.7 ~ 0.8 when x — > 1, al- 
lowing partly negative polarization of u quark, contrary 
to the result of 1 from the other three models. Then 
as for As v {x)/s v (x), SU(3)-extended BBS model and 
quark-diquark model give results around —0.5, against the 
positive value given by pQCD and brief pQCD. And as 
for s v (x)/u v (x), pQCD, brief pQCD and SU(3)-extended 
BBS models give a similarly positive value, while quark- 
diquark model predicts zero. 

Fig. 5 presents results of S~. The behaviors for 
As v (x) I s v (x) and Ad v (x) / d v (x) are very similar to those 
of Au v (x)/u v (x) and As v (x) / s v {x) in Fig. 3, respec- 
tively. And as for d v (x)/s v (x), it is shown that quark- 
diquark model get a slightly different behavior from that 
of the other models. 

The results of n, £~ and S° are not presented here, 
because they can be easily obtained from p, S + and S - 
respectively, according to u-d isospin symmetry. For ex- 
amples, u n {x) = dP{x) and dP (x) = vP (x). 

One can also use other phenomenological parameteriza- 
tions to study octet-baryon PDFs from SU(3) relations, 
and in Ref. [16] two sets of octet-baryon PDFs were ob- 
tained from polarized PDFs paramctcrizations of the pro- 
ton. 

Besides, our calculation of quark fragmentation func- 
tions also gives predictive results. According to our BBS- 
based model and the Gribov-Lipatov relation, we get 
DP{z)/D p d {z) -> 6.252 and D£(z)/D£(z) H- 1.648 when 
z — > 1. But in the statistical approach in Ref. [T5], the 
calculation gives DP(z)/D p d (z) -> 1 and L> A (z)/L» A (z) -> 
25.217. We can see that the difference is very large both 
for the proton and A. We can also compare with AKK 
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(Albino-Kniehl-Kramer) phenomenological parameteriza- 
tion |26| of the proton and A fragmentation functions, and 
the calculation results are both corresponding to the 
above two quantities. Thus further discriminations be- 
tween different predictions arc needed by experiments. 

Finally, we would like to estimate the effect introduced 
by SU(3) breaking and Mclosh-Wigncr rotation [37] to 
the SU(3)-extended BBS model. In the quark-diquark 
model, this effect comes from different masses of u, d and 
s quarks, and the correction is roughly estimated to be of 
about 10%. We therefore expect that such effect may bring 
about the same size correction in the SU(3)-extended BBS 
model. It means that our results are qualitatively reason- 
able, with some freedom to improve quantitatively when 
more data are available. 

Conclusion. — In this paper, we calculate parton 
distribution functions (PDFs) of octet baryons through 
SU(3) flavor symmetry, starting from PDFs of the pro- 
ton in the BBS statistical model. We make comparisons 
among different models, and find that our combination of 
the BBS statistical model and SU(3) flavor symmetry is an 
optional way to explore hadron structures. We also calcu- 
late the fragmentation function (FF) ratios, DP(z)/D^(z) 
and D£(z)/D£(z), at the z — > I region, based on the phe- 
nomenological Gribov-Lipatov relation. The results are 
compared with those from predictions of other models. 
We find our results arc reasonable and the method can be 
extended to the studies of baryon structures. 
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Table 1: PDFs of octet baryons in the quark-diquark model. The unit of m D (D — V or S) is MeV, and m q = 330 MeV for u, 
d quarks, m q — 480 MeV for s quark. 
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Table 2: Parameters for PDFs of octet baryons in pQCD. 
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Table 3: Parameters for PDFs of octet baryons in the brief pQCD. 
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Table 4: PDFs of octet baryons from SU(3) flavor symmetry. 
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Fig. 1: (a) Au(x)/u(x), (b) Ad(x)/d(x), (c) d(x)/u(x) for the proton. The full, dotted, dashed, and thick full lines are results 
for quark-diquark model, brief pQCD model, pQCD model, and BBS model, respectively. 
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Fig. 2: (a) Au(x)/u(x), (b) As(x) /s(x), (c) u(x)/s(x) for A. The full, dotted, dashed, and thick full lines are results for 
quark-diquark model, brief pQCD model, pQCD model, and SU(3)-extended BBS model, respectively. 
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Fig. 3: (a) Au(x)/u(x), (b) As(x)/ s(x), (c) s(x)/u(x) for E + . Models are identical to those in Fig. 2. 
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Fig. 4: (a) Au{x)/u{x), (b) As{x)/s{x), (c) s(x)/u(x) for E°. Models are identical to those in Fig. 2. 
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Fig. 5: (a) As(x)/s(x), (b) Ad(x)/d(x), (c) d(x)/s(x) for H~. Models are identical to those in Fig. 2. 



p-6 



